Sensitivity to thermal noise of atomic Einstein-Podolsky-Rosen entanglement 
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We examine the prospect of demonstrating Einstein-Podolsky-Rosen (EPR) entanglement for 
massive particles using spin-changing collisions in a spinor Bose-Einstein condensate. Such a demon- 
stration has recently been attempted by Gross et al. [Nature 480, 219 (2011)] using a condensate 
of *^Rb atoms trapped in an optical lattice potential. For the condensate initially prepared in the 
[F, mp) = (2, 0) hyperfine state, with no population in the niF = ±1 states, we predict a significant 
suppression of the product of inferred quadrature variances below the Heisenberg uncertainty limit, 
implying strong EPR entanglement. However, such EPR entanglement is lost when the collisions are 
initiated in the presence of a small (currently undetectable) thermal population Tith in the mp = ±1 
states. For condensates containing 150 to 200 atoms, we predict an upper bound of nth — 1 that 
can be tolerated in this experiment before EPR entanglement is lost. 

PACS numbers: 05.30.Jp, 03.75.Hh, 05.70.Ce 



Entanglement has proven to be "the characteristic trait 
of quantum mechanics" as first coined by Schrodinger [l| . 
It forms the foundations of quantum information theory 
and quantum computing. Further, in interferometry en- 
tanglement enables measurement precision to surpass the 
standard quantum limit [2|. This is particularly impor- 
tant in atom interferometry [3| as atom flux is gener- 
ally limited. However, the most important foundational 
trait of entanglement comes with its role in the Einstein- 
Podolsky-Rosen paradox (EPR) [J| , which pertains as to 
whether quantum mechanics can be considered a com- 
plete theory or not [5||. The EPR paradox for continuous- 
variable quadrature observables [6| (which are analogous 
to the position and momentum observables originally dis- 
cussed by EPR) has been demonstrated in optical para- 
metric down-conversion [7| and most recently attempts 
have been made to demonstrate the paradox with en- 
sembles of massive particles generated by spin-changing 
collisions in a spinor Bose-Einstein condensate (BEC) [8|. 

In this paper, we seek to provide a theoretical treat- 
ment of the recent experiment by Gross et al. [8| which 
reported entanglement, or quantum inseparability, of two 
atomic ensembles produced by spin-changing collisions 
in a *^Rb BEC. For the BEC initially prepared in the 
{F,mp) = (2,0) hyperfine state, the collisions produce 
correlated pairs of atoms in the mp = ±1 sublevels. The 
authors observed that the resulting state was insepara- 
ble, though a measurement of a stronger EPR entangle- 
ment criterion was inconclusive. A normalized product 
of inferred quadrature variances of 4 ± 17 was reported, 
whereas a demonstration of the EPRparadox requires 
this quantity to be less than unity |6|,|9|. 

The short-time dynamics of the spin-mixing process, 
for a vacuum initial state of the mp = ±1 atoms, is simi- 
lar to that of a spontaneous parametric down-conversion 
in the undepleted pump approximation. This paradig- 
matic nonlinear optical process is known to produce an 
EPR entangled twin-photon state that can seemingly vio- 



late the Heisenberg uncertainty relation for inferred opti- 
cal quadratures [6| . Such a violation has previously been 
observed in 1992 by Ou et al. [7|. Due to the inconclu- 
sive nature of an analogous measurement of matter- wave 
quadratures in Ref. J8|, we seek to perform a theoretical 
analysis of spin-changing dynamics and calculate vari- 
ous measures of entanglement in experimentally realistic 
regimes. In particular, we focus on the sensitivity of EPR 
entanglement to an initial population in the mp = ±1 
sublevels with thermal statistics. In the optical case this 
question is argued to be irrelevant as at optical frequen- 
cies and room temperatures the thermal population of 
the signal and idler modes is negligible, allowing us to 
safely approximate them as vacuum states. However, 
these considerations are inapplicable to ultracold atomic 
gases. Accordingly, when interpreting experimental re- 
sults care must be taken in differentiating spin-mixing 
dynamics initiated by vacuum noise from that initiated 
by thermal noise or a small coherent seed [10|. To this 
end, our modelling of the experiment of Gross et al. [8| is 
more consistent with a small thermal population in the 
mp ~ ±1 sublevels, rather than a vacuum initial state. 
From a broader perspective, our results for this widely 
applicable model highlight the generally fragile nature of 
atomic EPR entanglement to thermal noise, demonstrat- 
ing that future experiments must be refined to overcome 
this problem. 

The experiment of Ref. [i| starts with a BEC of ^"^Rb 
atoms prepared in the {F, mp) = (2, 0) state and trapped 
in a one-dimensional optical lattice. The lattice potential 
is sufficiently deep to prevent tunnelling between neigh- 
bouring wells and to also restrict spatial dynamics in each 
well to the ground state. Hence we may treat the conden- 
sate in each well according to the single-mode approxima- 
tion [111, in which the field operator ijji{r) for each com- 
ponent i = 77iF = 0, ±1, ±2 is expanded as ijji(r) = (j)i(r)di, 
where (f>i{r) is the ground state wavefunction and cii is 
the respective bosonic annihilation operator. We further 



assume (/'i(r) = (/)(r) for all components |ll| . 

A quadratic Zeeman shift and microwave dressing of 
the mp = state is employed to energetically restrict 
the spin-mixing dynamics to the mp =0, ±1 states [8|, 
and so for short time durations we may map the spin-2 
system to an effective spin-1 Hamiltonian [1^ of the form 
H = Hind + Hc\ + Hz , 

Hind = hg{alala-iai + a|a^iaoao), (1) 

Hd = fig{honi+non-i), Hz = -q{ni + h-i) , (2) 

where fii = ajai is the particle number operator and 
i ~ 0,±1 are referred to, respectively, as the pump 
and signal/idler modes from herein. The inelastic spin- 
changing collisions are described by -ffinci, and the re- 
maining elastic s-wave scattering terms are grouped in 
Hd, where g is the coupling constant associated with 
s-wave collisions 



12l . Il3fl . The quadratic Zeeman effect 



(Hz) due to the interaction with the magnetic field is 
parametrized via q, and we have ignored the linear Zee- 
man effect due to the conservation of the relative atom 
number between the signal and idler modes |13| . 

Simple analogies between the states of the signal and 
idler modes in spin-changing collisions and optical para- 
metric down-conversion consider only Hind in the unde- 
pleted pump approximation, however, competing mean- 
field {Hei) and Zeeman (Hz) effects lead to additional 
dynamics )14| due to dephasing. The full Heisenberg op- 
erator equations of motion are given by 



dr 
da±i 

dr 



2a_iaiaJ + {hi + n_i) ao 
agoL + {ho - q/g) a±i 



(3) 

(4) 



spinor c omp onents that decelerates the pair-production 
process [1J|. To prevent phase mismatch in the short- 



where we have introduced t = gt as dimensionless time. 
We see that the phase accrued in the a±i modes grows 
oc (no — q/g) whilst for the oq mode the phase grows 
ex {hi -\-h-i). In the short-time undepleted pump ap- 
proximation [l5|, this is equivalent to a phase rotation 
a±i — >■ a±iexp[z (A^o — q/g) t], where A^o = (fio(O)) is the 
initial population of the vip ~ component. This rota- 
tion leads to a dynamical phase mismatch between the 
: omp c 

. . s 

time limit one can choose q = gN^ in which case Eqs . 
©-([I]) reduce to those of resonant down-conversion |15J . 
We first analyze the spin-changing dynamics for the 
case of a vacuum initial state for the signal/idler modes, 
and a coherent state |ao(0)) for the pump mode with ini- 
tial number of atoms A^o = |ao(0)p. This case can be 
treated in a straightforward manner (see, e.g., Ref. [9|) 
by diagonalizing the full Hamiltonian in the truncated 
Fock-state basis and solving the Schrodinger equation 
(see also [16|). Figure [1] (a) shows the population dynam- 
ics of the signal and idler modes, for different initial atom 
numbers A^o and the quadratic Zeeman term tuned to the 
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Figure 1. (a) Fractional population nj^i{T)/No of the sig- 
nal/idler modes [where n±i(r) = {a[^^{r)a±i{T))\ as a func- 
tion of the dimensionless time r, for vacuum initial state and 
different initial number of atoms in the pump mode, Nq. The 
quadratic Zeeman term is phase-matched to q = gNo in all 
cases, except for the grey solid line which is shown for com- 
parison for q = and Ao = 175. The vertical dotted line indi- 
cates the measurement time r' = 0.0073 used in Ref. Q. (b) 
Same as in (a) but with thermally seeded populations in the 
signal/ idler modes (assumed to be equal to each other), for 
Ao = 175. The grey dashed lines show the analytic predictions 
in the undepleted pump approximation. 



phase-matching condition q = gNo. Setting q ~ elim- 
inates the Zeeman shift and we observe (grey solid line) 
significantly slowed dynamics due to phase mismatch. 
For reference, we also mark the experimental measure- 
ment time of Ref. |8|, r' = 0.0073, corresponding to the 
reported value of the squeezing parameter r = A'o''"' — 2 
[ig, evaluated for A^o = 175. 

We next analyze the case of an initial thermal seed in 
the signal/idler modes, with an equal average number of 
atoms nth in both modes. To simulate the dynamics in 
this case, we use the truncated Wigner method (Ref. |l7| 
gives simple prescriptions on how to model various ini- 
tial states in the Wigner representation). Figure [1] (b) 
illustrates that the presence of the thermal seed accel- 
erates population growth, however, it does not signif- 
icantly effect the maximal depletion of the BEC. The 
numerical results in Fig. [T] (b) are compared with the 
analytic predictions of the simple model of parametric 
down-conversion in the undepleted pump approximation, 
n±i{T) = sinh^(A^o ''')[! + 2nth] + hth [l3|- As expected, 
we find good agreement between the numerical and an- 
alytic results in the short-time limit. We also conclude 
that as far as the mode populations are concerned, the 
experimental measurement time t' = 0.0073 is not too far 
from the regime of validity of the simple analytic model, 
at least for nth ^2. This conclusion, however, cannot be 
carried through to other observables, such as entangle- 
ment measures analysed below. 

Central to this paper is an investigation into the pos- 
sible demonstration of the EPR paradox as outlined in 
Ref. [Si . In the context of continuous- variable entangle- 
ment, this is equivalent to the seeming violation of the 
Heisenberg uncertainty relation for inferred quadrature 



variances [a, |9| . In the normalised form this EPR entan- 
glement criterion can be written as 



T, 



{l-{a]a,)/m)y 



<1, 



(5) 



where the inferred quadrature variance for Xj (and sim- 
ilarly for Yj) is given by [6|| 



^inf-'^i 



{{AX.f) 






(6) 



with AXj = Xj — (Xj) and i,j — ±1. The general- 
ized quadrature operators here are defined as Xj (6) — 
ia]bje''^ + h]aje'~'''')/{h]bj)^/^ 0, where the operators h^ 
represent the local oscillator fields required for homodyne 
detection of the quadratures, and by convention we de- 
note Xj=Xj (0) and Yj=Xj (7r/2). 

In Fig. [21 (a) we show the results of calculation of 
a phase-optimized |13j EPR entanglement parameter T 
(with T_i = Ti = T due to the symmetry of the sig- 
nal/idler modes) for the signal/idler modes initially in 
a vacuum state. We see that strong EPR entangle- 
ment (T < 1) can be achieved for a large experimen- 
tal time frame, up to r ~ 0.01; more specifically, we 
predict suppression of the optimal EPR entanglement 
of at least 90% below unity for all relevant total atom 
numbers (ranging from 150 to 200) at t' = 0.0073. Un- 
like the simple undepleted pump model, which predicts 
T = cosh^ {2NaT) and hence indefinite suppression of 
the EPR criterion [15|, EPR entanglement in the full 
model is eventually lost due to a combination of back- 
conversion (1+1) + | — 1) -^ |0) + |0)) and the loss of co- 
herence in the pump mode. 

The significant error margin (T = 4 ± 17) of the ex- 
periment [8| notwithstanding, our results predict that a 
strong EPR violation should have been observed if the 
signal and idler modes were indeed generated from an 
initial vacuum state. In light of this result, we now dis- 
cuss the possible presence of stray or thermally excited 
atoms in the signal/idler modes and the effects such seed- 
ing can have on entanglement and particularly the EPR 
criterion. The results of calculation of the EPR entan- 
glement parameter T for an initial thermal seed of rlth 
in both modes are shown in Figs. |2|(b)-(d). We find the 
introduction of a thermal seed reduces the strong corre- 
lation between the signal and idler modes, leading to an 
eventual loss of EPR entanglement. For an initial num- 
ber of atoms in the pump mode ranging between 150 to 
200; EPR entanglement is lost already for fith ^ 1- Di- 
rect experimental detection of stray atoms at such a low 
population level is beyond the current resolution of ab- 
sorption imaging techniques [8[. More generally, our nu- 
merical results show that the maximum fith that can be 
tolerated while preserving the EPR entanglement scales 
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Figure 2. (a) Evolution of the EPR entanglement parameter 
T for the same situation as in Fig. [l] (a). The EPR crite- 
rion corresponds to T < 1 (dashed horizontal line). The inset 
shows the evolution of the optimal phase angle of the local 
oscillator 6o{t) for each A*'o [13l |. (b) Evolution of T for ther- 
mally seeded signal/idler modes and A^o = 175. The experi- 
mental measurement time r' = 0.0073 is shown in (a)-(b) as a 
vertical dotted line, (c) Time-optimized EPR parameter Tmin 
as a function of fith, for different A^o- The respective grey lines 
are the analytic predictions from the undepleted pump model. 
The grey line with squares shows Tmin for TVo = 175, but as- 
suming that the seeds are in a coherent state (sharing initially 
the same phase as the pump mode) with average populations 
of |a±i(0)p = nth- (d) Same as in (c), but as a function of 
A^o, for three different thermal seeds fith- 



|13| . For comparison, seeding the signal and idler modes 
with a coherent state |lO| of similar population does not 
have such a dramatic effect on EPR entanglement [see 
the grey line with squares in Fig. [2| (c)]. 

To further highlight the high sensitivity of EPR en- 
tanglement to initial thermal noise we contrast it with 
two other weaker measures of the nonclassicality of the 
state: two-mode quadrature squeezing and intermode en- 
tanglement in the sense of inseparability, which were the 
main focus of Ref. J8|. The two- mode quadrature vari- 
ances are defined as X± (6) = Xi (9) ± X-i (9), with 
A^X_(0) < 2 corresponding to two- mode squeezing |15J . 
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Figure 3. (a) Two-mode quadrature variances A^X± (6) at 
r' = 0.0073 as functions of the piiase angle 6 — do, for vacuum 
(solid lines) and thermally seeded (dashed lines) signal/idler 
modes; A^o = 175 in both cases, (b) Time-optimized minimum 
of A'^X-{6o) as a function of nth, for different A'^q. (c) Same 
as in (b), but as a function of A'^o, for different nth- 



i.e., suppression of fluctuations below the level dictated 
by a minimuni uncertainty state. We plot the results 
of our numerical calculations of quadrature variances in 
Fig. [3] (a) . From these results we observe that the mea- 
surements of Ref. [8| do not agree with the amplitude of 
the oscillation that we find for an initial vacuum state 
(solid lines). Rather they suggest the presence of a small 
thermal seed of nth — 1 (dashed lines). Further, calcula- 
tion of the minimum of A^X_ [Figs. [3](b)-(c)] highlights 
that two-mode squeezing is preserved for thermal seed 
populations up to nth — 1.7, which is consistent with our 
interpretation of the measurements reported in Ref. [8| . 
Next we define the sum of single-mode quadrature vari- 
ances as ^ A^ = 2(A^Xi -I- A^Yi) and the sum of two- 
mode quadrature variances, X) ^2 = A^X_ + A^y+. 
(Following the treatment of Ref. [8| we calculate the 
single-mode quadrature variances with canonical defini- 



tion of quadratures, Xj{9) = clje -f- 



ay 



. ) Insepara- 



bility of the produced mp ^ ±1 pair-entangled state is 
equivalent to X^^i/S^i < 1 (iSl- Figures S] (a) and 
(b) demonstrate that this measure of entanglement is far 
less sensitive to the presence of a thermal seed in com- 
parison to the stronger criterion of EPR entanglement. 

In conclusion, we have demonstrated that for an initial 
vacuum state in the signal/idler modes a strong suppres- 
sion of the EPR criterion can be achieved in the param- 
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Figure 4. (a) Time-optimized inseparability criterion for the 
quadrature entangled state, quantified via ^A^/X^^i < 1, 
as a function of nth, for different A'o. (b) Same as in (a), but 
as a function of A'^o, for different nth- 



eter regime of Ref. [8|, most importantly including the 
experimental measurement time of r = 0.0073. How- 
ever, we also establish that the strength of EPR entan- 
glement depends crucially on the nature of the initial 
spin-fluctuations. Specifically, we predict that for a pump 
mode of initially 150 to 200 atoms, a thermal initial seed 
of nth — 1 is sufficient to rule out EPR entanglement. 
Weaker measures of entanglement, such as inseparabil- 
ity, are still possible to observe as these are far more 
robust to thermal noise. Importantly, our results suggest 
that the measurement of EPR entanglement can serve as 
a sensitive probe of the initial state which triggers the 
pair pro duction process, beyond measures employed in 
Ref. [10| . This understanding of the sensitivity of EPR 
entanglement to initial thermal noise will hopefully lead 
to refining of spin-mixing experiments towards demon- 
stration of the EPR paradox with massive particles. We 
expect our findings to be also relevant to related pro- 
posals based on molecular dissociation [19|, condensate 
collisions 20l423|. and optomechanical systems [24|. 
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Model Hamiltonian. In Eqs. (1) and (2) of the 
main text, the coupling constant g is given by g = 
^(3^4 4- 452)/rf^r |0(i")|'*, where gp = Anfi^aF/m de- 
scribes s-wave scattering with total spin F^ characterised 
by scattering length ap [13 ■ For comparison, for an ac- 
tual spin-1 system, the coupling constant would be given 
by g= £^ / dh |(/)(r)|^, where gp = A-Kh^aF/m. 

In our representation of Hd we have used the fact 
that the relative number difference of the signal/idler 
modes [hi — n_i) is conserved due to the symmetry of 
the spin-changing collisions. The interaction with the 
magnetic field is described by Hi = —p{ni — n_i) — 
q {fii -|- 7i_i), where the linear and quadratic Zeeman ef- 
fects are parametrized, respectively, by p = ^bBo/2 and 
q = p'^/huups [251, with Ci;HFs/27r ~ 6.835 GHz being the 
hyperfine splitting frequency of ^^Rb [26| . Due to relative 
number conservation we may ignore the linear Zeeman ef- 
fect as a constant of motion which amounts to a constant 
phase shift and therefore does not affect the dynamics of 
spin-changing collisions. Furthermore, we may also ab- 
sorb the effects of the microwave level dressing (used by 
Gross et al. [g|) into a redefinition of the parameter q, as 
it simply leads to a shifting of the mp = energy level 
with respect to mp = ±1. 

Phase Optimization. Phase accrued due to Hq\ + 
Hz leads to a drifting in the phase relation be- 
tween the local oscillator and the signal/idler modes. 
This means that the conventional quadrature choice 
of XjiQ) and A'j(7r/2) may not measure the optimal 
violation of the EPR criterion. The accrued phase 
of Arg[(aJ(r)aJ(T)a_i(T)ai(T))] = 20^){t) can be can- 
celled by shifting the phase of the local oscillator clq — )• 
aoexp[i6'o(T)]. Hence the optimal choice of EPR quadra- 
tures become Xj(Oq{t)) and Xj(9q(t) + 7r/2). 



Generalized Quadrature Operators. Our choice 
of generalized quadrature operators [9| varies from the 
canonical form Xj{6) = aje^'^^ + aje'^ [l5| as it does 
not assume perfectly coherent, strong local oscillators. 
Instead, it takes into account the fact that in this ex- 
periment the local oscillators were derived, just before 
the measurement time, from a partially depleted and 
already incoherent pump mode by splitting it into two 
halves with a rf 7r/2 pulse [8|. This is equivalent to an 
optical beam splitter in which the output is given by 
bj = (ao + avac)/v^, where Ovac represents the vacuum 
entering the empty port of the beam splitter. As a con- 
sequence of this, Eq. (5) can be simplified by replacing 
the expectation value {b;bj) with 4(aJao). 

Undepleted Pump Approximation To invoke the 
undepleted pump approximation, we assume that the 



pump mode is initially in a coherent state with an ampli- 
tude ao(0) = \/l^ (which we choose to be real without 
loss of generality) and that it does not change with time. 
By additionally choosing the quadratic Zeeman effect to 
be phase- matched {q = gNo), we can reduce the model 
Hamiltonian to that of optical parametric-down conver- 
sion [15|, H = fi,x(ajaLi + h.c), in which x = ff-^o- The 
Heisenberg equations of motion following from this are 
da±i/dT - 



iNoal-i, where r = gt is a dimensionless 



time. Solutions to these equations are given by 



a±i(T) = cosh(iVoT)a±i(0) ~ ismh{NQT)aL{0), (SI) 



J 



which are physically valid in the short-time limit, gen- 
erally corresponding to less than 10% depletion of the 
pump mode occupation. 

Considering specific initial states for the signal and 
idler modes, these solutions can be used to calculate 
expectation values of various quantum mechanical op- 
erators and observables. For example, for a thermal 
initial state with an equal population in both modes, 
(ai(O)ai(O)) = (aLi(0)a_i(0)) = nth, the subsequent evo- 
lution of the mode populations is given by 

(aii(T)a±i(T)) = smh''{NoT)[l + 2nth] + "th, (S2) 

whereas the anomalous moments evolve according to 



(a±i(r)d;p(r)) = -i sinh(A^o''') cosh(iVoT)[l + 2nth\ 



(S3) 



Similarly, the EPR entanglement parameter is found to be given by 



(1 + 2nth)^ + W? [(1 + 2rith) cosh {2Not) - 1] [2 (1 + 2hth) cosh (2iVor) - 1] 



(1 + 2nth) cosh {2Not) -^[{1 + 2nth) cosh {2Not) ~ 1]' 



(S4) 



where we have assumed Nq ^ 1 . The minimum value of this quantity (with respect to time r) gives the maximal 
violation of the EPR criterion. 



^ min — 



x/2]Vo 



^ - (1 + 2?l.th) - 2^ (1 + 2nth)' - ((1 + 2nthf + l) V2m 



-2 



which is achieved at the optimal time 

1 



25^0 



-arccosh 



(S5) 



-i (1 + 2nth) + ^^J{l + 2nthf + 2No 



(S6) 



r 



From Eq. (jSSP we also determine the maximum allowable 
thermal population before EPR entanglement is lost. By 
numerical analysis we find a maximum seed of (nth) max — 
O.OSiVg^^ in the range 100 ^ iVo ^ 350. We find this 
compares reasonably with the results of full numerical 



simulations, which predict (nth)r 



OMN, 
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Furthermore we may also calculate the minimum two- 
mode quadrature variance, 

A^X^ = 2(1 + 2nth)[cosh(27VoT) - sinh(2AroT-)], (S7) 

and the inter-mode inseparability parameter (see main 
text). 



SA^/SA^ 



tanh(27Vor). 



(S8) 



Despite their limited applicability and the quantitative 
disagreement with the numerical results, the analytic pre- 
dictions of the undepleted pump approximation give use- 
ful insights into the qualitative aspects of different mea- 
sures of entanglement. For example, to leading order, 
Eqs. (|S4p and (|S7P predict, respectively, quadratic and 
linear growth of the EPR entanglement parameter and 
two- mode squeezing with the thermal seed fith- whereas 
the inter- mode inseparability, Eq. (jS8|) . is insensitive to 
nth- The predictions for EPR entanglement and two- 
mode squeezing are in qualitative agreement with the nu- 
merical results discussed in the main text, whilst we find 
weak linear growth with nth emerges for inter-mode in- 
separability due to depletion of the pump. These qualita- 
tive predictions highlight the lower tolerance and higher 
sensitivity of the EPR entanglement to thermal noise. 



